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Abstract 

We propose an objective and robust method to extract the electrical conductance of single 
molecules connected to metal electrodes from a set of measured conductance data. Our method 
roots in the physics of tunneling and is tested on octanedithiol using mechanically controllable 
break junctions. The single molecule conductance values can be deduced without the need for 
data selection. 

To determine the feasibility of devices based on single molecules and to assess their properties, a 
single or a few molecules have to be wired between two metal electrodes. This has become a reality 
only recently through different techniques such as scanning-probe microscopy, and mechanical and 
electromigration break junctions. 1-5 Using these techniques, the electrical conductances G of a 
broad range of molecular junctions have been determined 6-9 and gating of single molecules has 
been demonstrated. 10-13 
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These promising results are somewhat counterbalanced by the poor agreement in the conduc- 
tance values of single molecules reported by different groups. This disagreement reflects our present 
poor insight into the atomistics of single molecule junctions. To overcome junction-to-junction fluc- 
tuations, a statistical analysis has been proposed, in which G values of many junction realizations 
are represented in a histogram. This analysis has first been implemented in atomic junctions, 4 ' 14,15 
and has subsequently been used in metal-molecule junctions. 1,3,6 Peaks in the histogram point to 
preferred junction geometries. Evidence for the formation of few-molecules junctions was derived 
from the observation of a series of G values appearing at multiples of a fundamental single molecule 
value. The appearance of peaks in G-histograms is a very striking observation. However, to re- 
solve these, data selection schemes have been applied. 16-18 This situation is unsatisfactory, because 
there is at present no generally accepted objective selection criterion. We address this important 
question in this Letter. 

As a test case, we have chosen octanedithiol junctions. 1,3,17-21 We use a mechanically control- 
lable break junction (MCBJ) 22,23 setup and acquire many conductance traces in succession. We 
compare conductance histograms, which were generated with and without data selection. We show 
that the conductance value assigned to a single molecule is robust, and that data selections do not 
help to improve the results. The most convincing representation is found in a histogram of logG 
rather than G. 

The measurements were performed at room temperature, in a liquid environment. Figure^ 
shows a schematics of our MCBJ setup. 24 As flexible substrates, we use electrically isolated 
stainless steel sheets, over which gold leads are fabricated by electron beam lithography. The 
scanning-electron microscope (SEM) image in Figure^, shows the suspended 100 nm- wide region 
in the center of the gold leads. A flexible cell in the top of the substrate guarantees that the 
leads are always immersed in liquid, (see Figure ^) . The substrate is held by two supports at 
the periphery and a push-rod pressing from below. Bending the substrate results in stretching the 
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suspended Au bridge, which shrinks (b) until it breaks (c) and a gap of size d forms. The change 
in gap-size Ad is related to the vertical movement of the push-rod Az by an attenuation factor 
a = Ad/ Az — 1.6 — 4 x 10 -5 . 24,25 The push-rod is moved at a velocity v z = 30/im/s, so that the 
two Au leads separate at 0.5 — 1.2nm/s. 

We apply a constant bias voltage of 0.2 V, and record the variation of the current I through 
the junction during repeated open-close cycles. The current was measured with a current-voltage 
converter, which can automatically adjust its gain between 0.1 and 100V//iA. This allows us to 
register the conductance variation during the whole process, starting from the fused Au junction 
with G > Go := 2e 2 /h, until the formation of single molecule junctions, with conductance values 
orders of magnitude lower. 

The formation of a metal- molecule bridge evolves in several steps. First, the Au-bridge gets 
thinner (Figure^ 3 ), until a rather stable single-atom contact is established. A plateau in the G{z) 
curve is then expected around Go- When the atomic contact is finally lost, the conductance de- 
creases strongly. This decrease may be interrupted if a metal- molecule bridge is formed (Figure^). 
In that case, another plateau in the G(z) curve is anticipated. 3,26 Similar to atomic junctions, this 
metal-molecule-metal bridge holds via its chemical bonds the two sides together and postpones the 
breaking open of the Au electrodes (Figure ^i). 

To explore this process, we have performed groups of 100 consecutive open-close cycles for five 
different samples, both in pure mesitylene, and in a 1 mM solution of octanedithiol in mesitylene. 
In FigureEl we show representative G{z) curves during opening of the bridge. Whereas the curves 
in the main panel focus on values in the low conductance regime, i.e., at G ~ 10~ 4 Go, the inset 
shows data around G « Go, corresponding to the single-gold-atom contact. Figure [5] shows how 
the shape of the conductance curves G{z) is modified by the presence of octanedithiol molecules. 
Whereas G(z) decays in an exponential fashion in the pure solvent (curves to the left of the dotted 
line) , distinct plateau features may appear in octanedithiol containing solution. These plateaus are 
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the signature of the formation of single (few) molecule junctions. In some curves, jumps between 
plateaus at different G can be seen. In those cases, the molecular junction reorganizes, and the 
number of bridging molecules may change. In Figure [3] the first two curves (blue) to the right of 
the vertical dotted line, are rather "clean" . In contrast, the three next curves (red) are quite noisy 
just before plateau formation. This suggests that there is a large degree of molecular movement 
in the junction, until the octancdithiol molecules eventually lock between the leads. Finally, some 
G(z) curves measured in the presence of molecules do not display plateaus (last three curves in 
green). In this case, no stable single-molecule bridge has been formed. Such traces correspond to 
approximately 50% of the curves. 

Next, we focus on the statistics of our measured data. For each sample, we take all 100 con- 
ductance traces G(z), and determine the probability with which a particular G- value is measured, 
Ng(G). This is depicted in the conductance histograms of Figure^ (bin size: AG = 4 x 10 -6 Go). 
Whereas Na(G) decays smoothly in the pure solvent, there are distinct peaks appearing in the 
octanedithiol case (indicated by arrows). This suggests that particular molecular configurations 
form with a high probability. However, the peaks in Figured ar e masked by a strong background. 
One may therefore wonder, whether a particular data processing method could improve the sharp- 
ness of the peaks. In the literature, different procedures have already been used, but they have 
not carefully been compared with each other. In forming histograms, the proposed procedures 
consist of (a) disregarding G(z)-curves that do not present clear plateaus, 16,27 (b) only using data 
points that belong to plateaus, instead of taking the whole G(z) data, 20 (c) only using average 
values derived from the data points belonging to conductance plateaus, and weighting these by 
the plateau length; 17 and d) using conductance jumps. 18 In focusing on the plateau values, these 
methods do effectively eliminate a background. However, they can be subjective, as they involve 
decisions as how constant the signal has to be to define a plateau, or where the plateau exactly 
starts and ends. 
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We propose here an alternative method, which does not make use of any data selection. We 
take all data, and only subtract a background that is adapted to the physics of the problem. This 
method is as powerful as all the previous ones and, most importantly, it is fully objective. In 
proceeding, we note that the conductance must contain a tunneling contribution. 24 The tunneling 
conductance G is exponentially dependent on the gap distance d, i.e., G oc exp(— 2nd). Here, 



K = y/2rri(f)/h is the decay constant, <j) the apparent barrier height, and m the electron mass. 
Furthermore, d — a(z — zq), where a is the attenuation factor of the MCBJ, 24:25 and zq is defined 
as z(d = 0). Rewriting this, we find lnG = — 2naz + constant. It seems therefore much more 
appropriate to plot histograms of In G rather than of G. 

Making use of this expression, we can now calculate which is the expected tunneling contribution 
in the conductance histograms. If we denote with Ng, A^nd an d N z the respective probabilities 
of measuring a certain value of G, In G and z, we may write 

N G (G) dG = N lnG (lnG)d\nG = -N z (z) dz. (1) 

Here, N z {z) = R/v Zl where R is the data acquisition rate, and v z is the velocity of the vertical 
push-rod. In our case, both these quantities are constant: R = 500 points/s, v z — 30 /jm/s. 
Solving eq^for N\ nG yields 

N lnG (lnG) = — — . (2) 
2v z na 

Consequently, N\ nG is constant, whenever <j) and a are constants. Hence, in a lnG- or logG- 
histogram, tunneling shows up as a constant background which is easily subtracted. In Figure |3j), 
we show a log G-histogram built from the data in FigureOt (bin size: Alog(G/Go) = 5 x 10~ 3 ). 
A constant background is indeed present for G < 2 x 10~ 4 Go for the pure solvent (blue line), for 
which tunneling is the only expected contribution. In contrast, clear peaks appear in the presence of 
octanedithiol. The log G-histogram representation is very powerful for another reason: it presents 
a full overview of the data. At a glance, both the single-atom Au contact peaks (G ~ Go) and the 
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molecules signal (G < 10~ 3 Go) are seen. Between 10 _2 -10 -3 Go (depending on the sample) and 
Go, there is almost no weight in the histograms. This indicates that the Au atoms retract quickly 
immediately after breaking the gold atom bridge. 

The tunneling background, which is constant in a log G- histogram, is inversely proportional to 
G in a G-histogram, the latter being the representation the literature focused on so far. Solving 
cq ^for Ng yields 

As can be seen in Figure^, this expression perfectly matches the G-histogram of the pure solvent. 
The blue- line backgrounds of Figure |3i and b correspond to the same R/(2 v z na). 

We can use this property to subtract the tunneling background for the histograms on dithiol 
molecules. To this end, we fit eq |21 to our data from below. This background is shown in Fig- 
ure^ (black-dashed line). The same R/ {2v z na) gives the black-dashed constant background in 
Figure^. Subtracting it from the data yields a corrected histogram which is guided by the physics 
of tunneling. The result of this subtraction is shown in gray in Figure 0] (main panel and insets), 
for two different samples. Figure corresponds to the data of Figure 01 From this analysis we 
conclude that junctions with conductance values at multiples of 4.5 x 10 -5 Go are more favorably 
formed. This number is then assigned to the conductance of a single Au-octanedithiol-Au bridge, 
Gi. 

We will next compare our background subtraction method with other approaches based on 
curve selection. This comparison is shown in Figure 0] The blue-line histogram has been obtained 
by taking only curves in which plateaus are apparent (i.e. the blue and red curves in Figure [2J. 
In the red dashed histogram only the points within a plateau have been used. The latter data 
selection scheme is highlighted in black in Figure The selection was done manually and no other 
treatment was applied. 

On comparison of these three histograms, it is quite striking that all exhibit the same key 



features. There are two, sometimes even three conductance peaks at multiples of the same G- value 
(i.e., G\ « 4.5 x I0 _5 Go). Particularly interesting is that the gray and blue-line histograms in 
Figure 0] are almost identical. One can conclude from this that the G(z) curves without apparent 
plateaus can, on average, be described by a tunneling dependence. The effective barrier height in 
this case is somewhat smaller than that in the pure solvent. The third, red-dashed histogram, in 
which only plateau values were considered also yields similar peak positions, but appears to have 
an even stronger background subtracted. This is expected as in this histogram the noisy signals 
away from the plateaus (as shown in the red curves of Figure [2J have been removed. 

From the histograms, we find a single molecule conductance G\ w 4.5 x 10 -5 Go- In literature, 
values ranging from 1 to 25 x 10~ 5 Go have been reported. 1,3 ' 17 ~ 21 Our value lies very close to the 
one found by Wandlowski et al. 20 It is also close to that of Steigerwald et al. 28 for octanediamine 
in trichlorobenzene (2 — 6 x 10~ 5 Go), which was obtained without the need of any data selection. 
This similarity is particularly remarkable considering the different bonding group of the molecules. 
Tao et al. 17 reported two groups of peaks, at multiples of Gl = 5.2 x f -5 Go and multiples 
of Gh = 2.5 x 10 _4 Go. They attributed these to two distinct microscopic arrangements of the 
molecule-S-Au bonds. Whereas the first value agrees well with our findings, we do not observe any 
other peak at higher conductance values. This is especially made clear by the log G-histogram of 
Figure 03. The different solvent used in their work could be a possible cause for the formation of 
the second group of peaks. However, Tao et al. 17 observed peaks at the same conductance values 
in different solvents. Another notable difference between the two experiments is the speed at which 
the junctions are opened: 40nm/s in the work of Tao et al., 17 and I nm/s in our case. We speculate 
that the change in speed could lead to the detection of different microscopic conformations. Finally, 
a given microscopic arrangement could also be favored in our symmetric MCBJ, in comparison with 
the more asymmetric junctions formed in a scanning tunneling microscopy (STM) configuration. 

From the above discussion, it is clear that a detailed analysis of conductance histograms is 
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required to gain insight in the microscopic formation of single molecule junctions. In the analysis 
methods employed so far, a data selection process has been used. In contrast, we demonstrate that 
a simple background subtraction scheme suffices. It is as powerful as any data selection scheme 
and, in contrast to the latter, it is objective. We emphasize that the statistical analysis is most 
conveniently performed in a histogram in which log G is represented. The log G representation 
allows a simple background subtraction and provides an overview from the single atom contact to 
tunneling. Moreover, the single (few) molecule conductance values show up in a much more striking 
manner. In addition, we conclude that the features appearing in the conductance histograms 
obtained with break junctions (in MCBJ or STM configuration) are robust and can be realistically 
attributed to the molecular signature in these junctions. 
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Fi gure 1 ! (a) Schematics of the MCB J principle with liquid cell and a SEM image of one underetched Au junction. 
(b)-(d) Principle of the formation of a metal-molecule- metal bridge during the breaking process. Starting from the 
fused Au leads (b), a molecule can lock between the leads (c). Under further stretching, the Au leads are deformed, 
while the Au-octanedithiol-Au junction stays intact (d). 
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Figure 2: Variation of conductance during the breaking process of a junction in pure mesitylene (left of the 
vertical pointed line), and in a solution of octanedithiol in mesitylene (right). The curves are shifted in z for clarity. 
In the presence of octanedithiol, 50 % of the curves present plateaus. From these, some are very clean (the two 
first ones from left - blue), and others are noisier (the following three ones - red). The remaining 50 % (the last 
three ones - green) show an irregular decay without plateaus. The plateaus have been highlighted in black. Inset: 
Examples of plateaus close to 1 Go, corresponding to one-atom gold contacts. 
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Figure 3: (a) Conductance histograms built from approximately 100 G(z) curves (Figure in pure mesitylene 
(red), and in a solution ImM of octanedithiol (grey). The arrows indicate the conductance peaks that appear when 
octanedithiol is added in solution. The blue and black-dashed lines show the best fit from below using a expression 
oc 1/G to both histograms. (b) Histograms of logG built from the same data as in (a). The blue and black-dashed 
lines correspond to the same R/ (2 v z na) values as in (a). Note that in the latter representation both atomic gold 
peaks and molecular peaks are observed. 
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Figure 4: Histograms built from a given group of G(z) curves, (a) and (b) show results for two different samples. 
The grey histograms were obtained considering all the measured conductance traces, and subtracting later the 
tunneling background (also shown in the insets) . The blue-line histograms were made with the curves that display 
clear plateaus (blue and red curves in Figure |5J. Finally, the red dashed histograms were built only with the G 
values which belong to plateaus (marked in thick black in Figure |5J. The Gaussian curves highlight the position of 
the peaks. 
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